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Abstract. We study the short-time existence and uniqueness of solutions to a coupled system of 
partial differential equations arising in mean field game theory. It has the generic form 

J —dtu — Au + H{t,x,m,Vu) = f{t,x,m) 

\ dtm — Am — div {mVpH{t,x,m,Vu)) — 0 

plus initial-final and boundary conditions. The novelty of the problem is that the Hamiltonian 
H{t,x,m,p) may take such forms as m~°‘\pY for some a > 0 and r > 1. Our main result is the 
existence of weak solutions for small times T so long as r is not too large, and uniqueness under 
additional constraints. The main ingredient in the proof is an a priori estimate on solutions to 
the Fokker-Planck equation. We also briefly consider existence and uniqueness of solutions to an 
optimal control problem related to mean field games. 

Keywords: mean field games, Hamilton-Jacobi, Fokker-Planck, coupled systems, optimal control, 
nonlinear partial differential equations 
AMS Classification: 35K61 


1. Introduction 


Let O C be a bounded domain, let T > 0 be given, and set Q := O x (0,T). Our purpose is to 
study 


( 1 . 1 ) 


(i) —dtu — Au + H{t,x,m,'Vu)=f{t,x,m) in Q 
< {ii) dtm — Am — div {m'VpH{t,x,m,Vu)) = 0 in Q 
(in) u{T,x) = g{x,m{T,x)), m(0, x) = mo(x) in 0 


equipped with Dirichlet, Neumann, or periodic boundary conditions. Here we assume that H{t,x, m, p) 
has a structure of which the canonical example is 


( 1 . 2 ) 


H{t, X, m,p) = 


\P\' 


rm^ 


■, r > 1, a > 0. 


The main result of this paper is the existence of short-time solutions for system 0 under the 
condition that r is not “too big” (see Section EH). 


System represents a mean field game with a congestion term. Heuristically, the function u 

denotes the value function of a representative player whose objective is to minimize 

(1.3) 


E 


L{s,X{s),m{s,X{s)),v{s,X{s))) + f{s,X{s),m{s,X{s)))ds+g{X{T),m{T,X{T))) 


Ut 


Research supported by NSF grant DMS-1303775. The author wishes to thank Alain Bensoussan for his mentoring 
and support. 
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over all trajectories X = X{t) of the controlled stochastic differential equation 
(1.4) dX{s) = v{s,X{s))ds + y/2dB{s), X{t) = x. 


The Lagrangian L is defined by taking the Legendre transform of H in the last coordinate, the 
canonical example being given by 


L(t, X, m, v) 


m \v\ 
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Here m{t, x) is the density of participants in the game, while a can be taken as a congestion 
parameter. Increasing a will increase the relative cost incurred by the representative player when 
moving at high speeds in areas of high population density. 


The main result of this paper is the existence and uniqueness of weak solutions of (ll.lji under 
general assumptions on the data. Our main contribution is to show that H{t,x,m,Xu) may have 
a highly singular dependence on m, provided that it does not blow up too quickly as Vtt becomes 
large. This condition is made precise in Section [2.11 below. 


To put this into perspective, we note that the vast majority of existence results for mean field 
games address only the case where H does not depend on m at all, that is, H = H{t,x,Xu). 
Using a priori estimates and compactness criteria, one can construct weak solutions using known 
results for quasilinear parabolic equations [nun]. The other option in this case is to view (HTD 
as a condition of optimality for the control of the Fokker-Planck equation (or, by duality, of the 
Hamilton-Jacobi equation). This was pointed out in [12] and exploited in jS] UHl EJ [5] to obtain 
new existence results even for systems with degenerate diffusion (that is, replacing the Laplacian 
— A with an elliptic differential operator which is not strictly positive). In such works, the loss of 
regularity is made up for by the fact that solutions can be constructed by passing to the limit on 
optimizing sequences, which, with the help of a priori estimates on Hamilton-Jacobi equations, can 
be shown to be compact in appropriately chosen functions spaces. See |5] for details. 

Up until now, very few authors have studied PDE systems of mean field games which did not fit 
this paradigm. The motivation for the present work is a result of Gomes and Mitake in |8|, where 
they prove the existence of classical solutions for a stationary version of (11.11) with the Hamiltonian 
given by the canonical example (II.2p with r = 2. Their proof relied on a priori estimates which 
emerge from the specific structure they consider. In particular, integration by parts yields a special 
cancellation which allows one to conclude that 1/m is bounded, after which everything else follows 
from classical methods for parabolic equations. 

A system which resembles the mean field game dnD where H has the form ([O]) appears in [4]. 
However, in that reference the problem is one of optimal control, and the PDE system studied 
there is a condition of optimality and not a game. Such optimal control problems arise in transport 
theory [7| and also have applications to traffic flow and congestion [2]. It is interesting that the 
PDE systems for such control problems resemble (ffTD. However, the fact that our case is not a 
control problem creates added difficulties which we must surmount. See [T] for related remarks on 
the comparison between mean field games and mean field optimal control problems. We will briefly 
visit the mean field type control problem corresponding to dni) in Section m 

Rather than considering only the canonical Hamiltonian given by (|1.2p . our main results cover 
much more general Hamiltonians which may depend on m in a highly singular way. The difficulty 
we face when constructing a solution is showing compactness in the Hamiltonian term. Classical 
methods for quasilinear parabolic equations (see m and references therein) do not apply here, 
since H{t,x,m,Xu) could blow up due to m rather than Xu. In this work we deal with this by 
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finding conditions under which a solution to o will satisfy the a priori estimate 

||"l||oo + ||l/"2.||oo < C, 

which then allows the classical methods to go through. Such an estimate can be proved for solutions 
of the Fokker-Planck equation provided that the vector-field generating the flow is sufficiently 
integrable. This is the main reason for assuming that H does not grow “too fast” with respect to 
the variable Vn : with this requirement, the energy estimates for mean field games yield that VpH 
is sufficiently integrable. 

Another work by Gomes and Voskanyan [9] studies a model similar to (ll.ljl . There they consider the 
existence of smooth solutions, for Hamiltonians of the form (11.21) having any subquadratic growth 
in the gradient variable but for a sufficiently small exponent on the density: that is, any 1 < r < 2 
by for a sufficiently small. The result holds, as in the present work, for sufficiently small times 
T. These two works can be seen as complementary insofar as (i) the present work does not study 
bounded, smooth solutions but rather weak solutions, and (ii) the present work studies existence 
for a in any range (but with restrictions on r). 

The author wishes to thank Diogo Gomes for helpful comments made during the preparation of 
this article. 

To conclude this introduction, we outline the structure of the paper. The section immediately 
following this introduction gives some notation, defines weak solutions and presents our main 
results, Theorems l2.2l and l2.4l Section[3]is devoted to proving the existence of solutions. The largest 
part of this section is devoted to obtaining a priori estimates on the Fokker-Planck equation-see 
Proposition 13.11 Section 0] is devoted to proving uniqueness of solutions under certain conditions. 
Finally, in Section [5] we briefly examine the problem of mean field type control. It is shown in 
particular that uniqueness of solutions for the control problem holds under much more general 
conditions than for mean field games. The assumptions in play in the present work make it easier 
to construct solutions of the mean field type control problem; we can expect that in future work, 
much more general results will be proved using different variational techniques. 

2. Preliminaries 

2.1. Notation and assnmptions. We denote by LP{Q) = LP((0,T) x Q) or just the Lebesgue 
space of p-integrable functions on Q, with norm || • \\lp or simply || • ||p. We denote by V the 
d-dimensional spatial gradient. Sobolev spaces will be labeled or when p = 2, with 

corresponding norms || • We will write || • for the norm in LP(0,T; T'^(H)), or more 

generally || • norm in VF*’^(0, T; The space of continuous functions on Q 

will be denoted C{Q), while C^{Q) denotes the space of all functions with s continuous derivatives. 
Analogous notation applies for function spaces on H. If X is any function space, then the space 
of all non-negative functions in X will be denoted Throughout this work the symbol C will 
denote a constant depending only on the parameters given by the theorem or proposition being 
proved; its precise value may change from line to line. 

We assume the following hypotheses on the data: 

(1) (Initial condition) mo S L°°(H), mo > c > 0, and /g mo = 1 

(2) (Goupling) f{t,x,m) and g{x,m) are measurable in {t,x), continuous and nondecreasing in 
m, and bounded below. Moreover, 

(2.1) fp^[t,x):= sup /(f,x,m), c/l(x) := sup gix,m) 

mg[0,L] 
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are both integrable for all L > 0. 


(3) (Hamiltonian) We assume that H = H{t,x,m,p) is measurable in {t,x), continuous in m, 
and both convex and continuously differentiable in p. Moreover, we assume that for some 
constants C > 0, A > 0, 


( 2 . 2 ) 

(2.3) 

(2.4) 


i^r 


H{t,x,m,p) > - . , 

C[m^ + m~^) 

X, m,p)\ < C{m^ + + |pr~^) 

VpH{t,x,m,p) ■ p — rH{t,x,m,p) > —C 


The exponent A appearing in four different places in (I2.2p - ()2.3p could in principle be four 
different exponents, but it is sufficient for our purposes to take A to be their maximum. 
The main point is that H and VpH can blow up polynomially in m or m~^, that is for m 
large or small. 


(4) (Growth of Hamiltonian) The exponent r setting the growth rate of H satisfies the condition 


(2.5) 


r < 


d + 2 
d+ 1' 


The canonical example presented in the introduction is H{t,x,m,p) = which satisfies 

the above hypotheses with X = a. However, we can allow much more general Hamiltonians, for 
instance those having the form 

H{t,x,m,p) = h{t,x,m)\p\^ 
where h is continuous in m and satisfies 

e{m^ + < h{t, x, m) < C{m^ + m~^) 

for some positive constants C, e, and A. Note that A may be large, so that h may be very large or 
quickly vanish when m is either large or small. 


2.2. Definition of solutions. We now define solutions of system m- 

Definition 2.1. A couple {u,m) G L^{Q) x L^(Q)_|_ is a weak solution of (ini) provided 

(i) m £ C{[0,T]] L^{Q)) and g{-,m{T)) £ 

(ii) (1 + m)| VpH(-, •, m, Vu)p, /(•,•, m) G L^(Q); 

(iii) the equations hold in the sense of distributions, 

(2.6) [ u{4>t-A(f)+ f H{t,x,m,Vu)(f = [ f{t,x,m)(f+ [ g{x,m{T))(f){T) 

J Q 'J Q Q J ^ 

for every (f £ C'“((0,T] x H) satisfying (j2.8h . and 

(2.7) / m(—— A(/)) + / mVpH{t,x,m,Vu) -Vf) = / mo(j){0) 

J Q J Q J 

for every cf £ C“([0,T) x H) satisfying (|2.8p . where the boundary conditions may be one of 
the following: 

periodic: cj) is — periodic in x, 

(2.8) Dirichlet: <f = Q on Ti := (0,T) x dXl, or 

Neumann: ^ = 0 on S := (0,T) x dQ. 
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A special remark is in order concerning the integrability condition [l + rn)\'S/pH{m,'S/u)'\^ G L}-{Q). 
The Fokker-Planck equation in ( 11 . 11 ) would still have a meaning in the sense of distributions given 
the weaker condition m\7pH[m,'Vu) G L^{Q). However, the stronger condition allows us to prove 
that weak solutions enjoy more regularity than given by the definition. See Proposition 14.21 Such 
extra regularity makes possible our proof of uniqueness given in Section HI 

We should also note that |VpFf(m, Vu)p G L^{Q) does not follow from H{m,'Vu) G L^{Q) by the 
assumptions on H. This is in contrast with the case where H does not depend on m. Indeed, if 
we removed the m dependence from (j2.3p . then H{m,\7u) G L^{Q) would imply \7pH{m,'Vu) G 
L’'/(’'-i)((5) C L‘^{Q). In actual fact, (12.3p does not allow us to deduce VpH{m,'Vu) G L‘^{Q) unless 
we already know that m and m~^ are bounded. 


2.3. Main results. The following two theorems constitute the main results of this paper. 

Theorem 2.2 (Existence). Assume the hypotheses given in Section \2.1l There exists a constant 
To > 0 such that ifT< Tq, then there exists a weak solution {u,m) of (ini), defined in Section [gJl 
below, such that 

(1) u,m ^ ^^([0, T]; L^(H)) with u{T) = g{m{T)) and m{0) = mo, m > 0; 

(2) m £ L‘^{0,T; H\n)) n L°°{Q), 1/m G and dm £ 

(3) u is bounded below. 

Furthermore, suppose that and gi, defined in \2.1l are bounded for all T > 0. Then u is bounded 
and locally Holder continuous in Q. 


In order to formulate our uniqueness result, we will need an additional hypothesis on the structure 
of the Hamiltonian. 

Assumption 2.3. ITe assume, in addition to the assumptions in Section \2.1\. that H is C‘^ in p 
and that for all m, m > 0 and p,p £ 

(2.9) — dmH{t, X, m,p)rh? + mVpH{t, x, m,p){p,p) + mrhp ■ dmVpH{t, x, m,p) > 0. 

A more concrete structure satisfying this assumption is given by 

(2.10) H{t,x,m,p) = h{t,x,m)\pf 

for some function h{t,x,m) which is measurable in {t,x), continuously differentiable in m, and 
satisfies 

4fr — 1) 

(2.11) 0 < —dmh{t,x,m)m < - -h{t,x,m) Mm > 0. 

r 

Note that if we take the canonical model from the introduction, then (|2.11l) implies 

h{t,x,m) = m~°‘ 0 < a < 4(r — l)/r. 

Theorem 2.4 (Uniqueness). In addition to the hypotheses given in Section fUTl suppose that H 
satisfies Assumvtion \2.3l Then the weak solution provided by Theorem \2.2\ is unique. 
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3. Existence 


The proof of Theorem 12.21 is based on a priori estimates and compactness. This section is divided 
into three parts, the first two of which address a priori estimates. First, in Proposition 13.11 we 
prove that solutions of the Fokker-Planck equation are a priori bounded both from above and away 
from zero, given a sufficient integrability condition on the vector field. The other a priori estimates 
are standard for all mean field games. Combined with Proposition 13.11 and known compactness 
results, we prove the existence of solutions to (ll.ljl . 


(3.1) 


3.1. A priori estimates on the Fokker-Planck equation. Consider the equation 

dtm — Am — V • (bm) = 0 in Q, 
m(0) = mo in 11 

with one of the following boundary conditions: 

periodic: fl = b, m are — periodic in x, 

(3.2) Dirichlet: m = 0 on S := (0,T) x dQ, or 

Neumann: ^ + mb • i/ = 0 on S := (0, T) x dQ 

where u is the outward unit normal vector on dQ. If 6 S L?‘{Q), we say that m G L‘^{0,T] (H)) n 

L°°((5)+ is a solution of (j3.ip - (l3.2l) provided that 


/ —mdt4> + Vm • V(/> + mb ■ S/cj) = / mo(j){0) 
Jo Jn 


corresponding to (13.2p . 


'Q 

for all test functions (p satisfying the boundary conditions from 

Proposition 3.1. Let mo G (17(11) be such that mo > 0, and let b G Lp‘{Q). Suppose m G 
L^(0, T; Ff^(ll))nT°°((5)_i_ is a solution of (I3.ip - ()3.2h . Suppose further that there exist non-negative 
measurable functions bi,... ,bN and real constants Pi,..., I3i\f such that 


N 


(3.3) 


|6| < m^'^bk G A: = 1, 


,iV. 


i=l 


Assume also that (ESP is satisfied. Then there exist constants To > 0, 

C(||"io||oo,{||m^''^fc||r/(r-i)}f=i) and C{\\mQ^\\oo, {\\m^’‘h\\r/{r-i)}k=i) 


such that ifT < Tq, 

then 

(3.4) 

sup 

te[o,r] 

(3.5) 

sup 

te[o,r] 


1 


m{t) 


< C(ll"ioioo,{||m^''6fc||r/(r-l)}f=l) 


Moreover, from dsap and (1331) we deduce 

(3-6) l|5t"i||L2(0,T;H-l(O)) + ll"^llL2(0,r;i^l(O)) < C'(||mo||oo, { ||r/(r-l)}^l)• 

Proof. First we will prove (13.5p . Assume that m is strictly positive; we will remove this assumption 
at the end of the proof. Multiply (13.11) by m~'^~^ and integrate by parts, then use Young’s inequality 
to estimate 


(3.7) f m{t) f f m < f 

J Q, J 0 J fl Jq 


m-o^ + 


q{q + l) 


f f m '^|6p. 

Jo Jn 


2 


2 









By the Sobolev embedding theorem we have 
i-T pT 


(3.8) ^ \\m <C ||V(m + \\m 

^ [ f dxdt + C f ||m“^(t)||^dt. 

4 Jo Jn Jo 


= C- 


With the previous estimate this implies 
(3.9) < C\\mQ^\\l^ + Cq^ 


m '^| 6 p 


/o Jo. 


N 'jp 

<C\\mQ^\\l^ + CNq^^ ! [ m-‘^bj 

“t Jo Jfi 


i^^Jo Jn 


N 


< C\\mQ^\\l^ + CNq‘^Y,\\md^bi\\l/^^_^)\\m ^^^(^+ 7 ,) 


i=l 


where r* := r/{2 — r), 7 * := max{2/3j(r — l)/r,0}, and where C is some constant depending only 
on the domain O and the final time T. 


\\q{l + 2/d) < ||"l IIl7>L| 11”^ \\jqjdq/(d-2)^ 


T g T dq/{d —2) ' 


Noting that 

(3.10) 
we can deduce 

(3.11) \\'m~^\\q{l+2/d) < ^Sl^/''g^/''max{l, ||mo^||oo, ||r*(g+ 7 )}^’^'^'^^ 

where 7 := max{ 7 i,... , 7 Ar} and K depends only on {||m^'' 6 fc||r/(r-i)}^i and constant parameters. 


-ii 


In order to use the Moser iteration method, we need to show that for some q large enough, \\m 
can be bounded by a constant depending only on the parameters as do C and K. Let us observe 
that (13.7|) can be used to deduce 


(3.12) 


f m ‘^{t)d. 

Jn 


lx + \\m ^(s 




ll^dq/(d- 2 )a* — 


ds < C + f f m dxds, 

Jo Jn 


where again C depends only on the parameters. We wish to bound the right-hand side. Choose 
A G (0,1) such that 

qX + {l- \)dq/{d-2) = {q + 7 )r*, 

that is, 


^ ^ dq - {d-2)r*{q + 'y) 


We also set 




2q 

d{r* - 1 - 1 - 'yr*jq) 


Note that /i —>• |(r* — 1) < 1 as g ^ 00 , hence for q large enough we have ^ < 1 as well. So 
estimate 


we 


(3.13) f m '"*^‘^~^^'^{t)dx < ( [ m ‘^{t)dx 
Ju \Ju 


m 


— dq! {d—2) 


(t)dx 


l-A 


m‘i{t)dx'^ \\m ^(t)||rfg/(rf_ 2 ) < C* 




+ \\m 




dq/(d—2) 
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where in this line C depends only on /i. Plugging this into (|3.12l) . we see that 

rt / f \ A/(i-/^) 

(3.14) 


f m ‘^{t)dx<C + C [ ( [ rn'^{s)da 

Jn Jo \Jn 


ds. 


Set /3 = A/(l — n) and H{t) = f* (f^m'^(s)dx)^'^^^ ds. This implies the ordinary differential 
inequality 

H{t) < C{l + H{t)f, H{0) = 0. 

In particular, we have 

H{t) < z - zzT— - -rn—r V t < := 


1 - c{u - 


C{(3 - 1) 


Hence if T is small enough, we have that H{t) and thus also \\m ^{t)\\q is bounded by a constant 
for all t G [0, T]. This will allow us to start the Moser iteration. 


We now dehne sequences {g„} and {zn} as follows. Set s := {1 + 2/d)/r*; by the assumption 
r < 1 + (d + 1)“^ we have s > 1. Fix some 

(3.15) 1 < u < s and qo > - 

s — u 

where go is large enough so that ||r*(ijo+ 7 ) ^ C*. Now set 

(3.16) Qn+i = QnS - 7, Zn = max{l, ||mQ ^lloo, V n > 0. 

Note that {qn} grows exponentially: 


(3.17) 

Apply (13.lip to get 

(3.18) 

Using the bounds on g„ this implies 


qou^ <qn<qos^- 

Zn+i < 


(3.19) log(z„+i) < 
where 


log(i7gg s 

gou" 


2 Xln\ 


+ 1 + 


7 




log(^;n) < 


A{n + 1) 


+ (1 + ^1 


4 _ niax{log(A:gg),log(s^)} ^ ^ 1_ 


By induction, one deduces from (j3.19p that 


qo 


qo 


(3.20) 


log(^n) il + ^) max{log(^o),l}, 


\k=0 


k=0 


which converges because r > 1. Setting Zoo '■= lim sup^^j^ we can conclude that 

(3.21) ||m“^||oo < Zoo <+00 

where Zoo is a constant depending only on 11, T, ||mQ ^||oo and \\'m^^bk\\r/(r-i)}k=i- Appealing once 
more to ()3.7p we get more precisely the estimate (13.5p . 


In the proof above, in order to remove that assumption that m is strictly positive, simply replace 
m in with m + e for e > 0 arbitrary. Then let e —^ 0 to get the result. 

To prove (13.4p . we use an analogous argument. In particular, multiplying (13.11) by we arrive 

at 

(3.22) - [ m{ty + -— ^ < - / m'l + f f m^\b\'^. 

q Jn Jo Jn q Jn 2 jg Jq 
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By the Sobolev embedding theorem and some calculations (cf. (13.91) 1 we obtain 


N 




i=l 


where in this case ji = max{—2/3j(r — l)/r, 0}. This leads, using the same argument as above, to 

(3.24) \\m\\q(^i+ 2 /d) < raax{l, ||mo||oo, \\'m\\r*{q+^)V^'^^‘^ 

where as before 7 = max{ 7 i ,... ,'yiy}. Then a similar iteration procedure as before leads to (13.41) . 


Finally, we deduce (j3.6p from (j3.4p and (13.51) as follows. Plug in the estimate (13.4p into p3.22p with 
g = 2 to see that m is bounded in T^(0, T ; Then observe that p3.3l) together with (13.4p and 

(j3.5l) imply that mb G We deduce from the Fokker-Planck equation that dtm is bounded 

in L 2 ( 0 ,T;i?~i(f 2 )). □ 


3.2. More a priori estimates. Here we collect some standard energy estimates for the mean field 
game system. 

Lemma 3.2 (Main energy estimate). Let f = f{t,x,s),g = g{x,s), and H = H(t,x,s,^) be given, 
and suppose {u, m) is a bounded solution of dni). Assume f,g, and H are bounded below by some 
constant —C. Then, suppressing the dependence on {t,x), we have 


(3.25) 


f g{m{T))m{T)+ f f /(m)m + (Vtt-Vpi/(m, Vu) —i/(m, Vtt))m + 
Jn Jo Jn 


10 Jn 


H{m, Vtt) 


< CII^Trolloo { / 5211)7101100+ / / /2||mo||oo + 

In Jo Jn 


where / 2 ||mo||oo 52||mo||oo defined in dH]). 

In particular, if mQ G L°° and / 2 ||moi|oo 52||mo||oo integrable, then 

(3.26) 

||/(m)m||ii^ + || 5 (m(r))m(r)||ii+||m(Vu-VpH'(m,V'u)-H(m,Vu))||ii^ + ||F(m,Vu)||ii^^ < C. 


Proof. We rely on the rather simple but important observation that 

(3.27) /("i) < /l + Yf{m)m + Cm, g{m) < gi + —g{m)m + Cm 

1j 1j 

We will apply this below with L = 2||mo||oo- 
Start with the energy identity. 


(3.28) f g{m{T))m{T) + f f f{m)m + {Vu-VpH{m,Vu) — H{m,Vu))m= ( u{0)mo, 

Jn Jo Jn Jn 

which is obtained by adding m times the Hamilton-Jacobi equation to u times the Fokker-Planck 
equation and then integrating by parts. Then using the lower bound on u, 

(3.29) 


f u{0)mo< [ (tt(0)-F C')mo < llmolloo [ ("(^(O)C) < ||mo||o 

g{m{T))+ [ [ f{m)- [ [ H{m,Vu) 

1 Jo Jn Jo Jn 


I u{T) -I I dtu + C 
n Jo Jn 

T 

+ C 


9 




Use (13.271) and plug into (|3.28p to get (I3.25p . 


As for (j3.26p . it suffices to observe that each of the terms in (|3.25p is non-negative: f{m)m and 
g(m{T))m{T) because / and g are monotone, Vu • Vp77 — H because H is convex in p, and H by 
hypothesis (12.21) . □ 

Corollary 3.3. Supose {u, m) is a smooth solution of ([HD. The following estimates hold for a 
small enough time T and for a eonstant C depending only on the data. 

(1) Estimates (j3.4p . (13.51) . and (13.61) hold. 

(2) \\VpH{m,Vu)\\r/(r-i) <C. 

(3) \\f{m)hl^ + Mm{T))U.<C. 

(4) \\dtu + Au\\Li^ < C. 

Proof. 1. As a result of (|2.3p . (12.21) . and then (|2.4I) . we have 

(3.30) \VpH{t,x,m,p)f/^^-^') < C{m^ + + \pf) 

< C{m^ + -|- H{t, X, m,p)) 

< C{m^ + + VpH{t, x,m,p) ■ p — H{t, x, m,p)). 

Let b = VpH{t, X, m, Vtt) and n = (1 + VpH{t, x,m,p) ■ p — H{t, x, Then we have 

|6| < + Cm^\ + Cm-^\ =: 6i + 62 + 63 + h^. 

Moreover, from (|3.26l) and the definition of v it follows that ||i < C. Therefore we can 

choose some jdi, i = 1,2, 3,4, such that < C. Hence we reach the desired conclusion 

by Proposition 13.11 

2. Combining the estimate ||m||oo + ||?tt-~^||oo < C from the previous step with ()3.30l) and (I3.26p . 
we conclude that ||Vpi7(m, Vu)||r/(r_i) < C. 

3. This follows from (I3.26P together with (|3.27l) . 

4. This follows by using the estimate ||i7(m, Vu)||i < C from (I3.26P together with part 4 in the 

Hamilton-Jacobi equation. □ 

Remark 3.4. ITe observe that hypotheses ()2.2p . (l2.3p . and (12.41) were not used directly to prove the 
a priori estimates in Corollary lg.3i but only to derive the estimate (I3.30p . This will be useful in 
the following section, where we show that a particular approximation of H which is Lipschitz in p 
satisfies ()3.30p . even though it does not satisfy (I2.2p . ()2.3p . and (12. 4p . 

3.3. Proof of Theorem 12.21 In order to construct solutions of dm), we first derive smooth 
solutions to an approximation of dm, and by compactness we show that they converge to weak 
solutions of dm- For the hrst step, we have 

Lemma 3.5. Let H{t,x,m,p) satisfy the following: 

\H{t,x,m,0)\ -(- \VpH{t,x,m,p)\ < C. 

Suppose f,g are also bounded. Then there exists a bounded (smooth in (0,T)j solution {u,m) to 

dm- 
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Proof. Set 


XL = {me Ci[0,T]-,L^{n))nL°^{Q) : \\m\\oo < L}. 

For any fj, G Xl, define G L^(0, T; n L°°{Q) to be the (unique) bounded solution of 


(3.31) 


—dtu — Au + H(t, X, m,Vu) = f{t, x, fi) 
u{T,x) = g{x,n{T,x)) 


We then set m = d*(^) to be the solution of 
(3.32) 


dtm — Am — div {mVpH (t, x, m, Vu^)) = 0 
m(0, x) = mo{x) 

By standard results on parabolic equations there exists a constant L depending on ||VpFf||oo, ||mo||cx) 
such that m is bounded by L uniformly in /r. For this L, it follows that is an invariant convex 
subset of (^([OjT]; L^(fl)) under the map <h, which is continuous and compact. By Schauder’s fixed 
point theorem, we conclude. □ 

We now define an approximation of system (11.11) satisfying the hypotheses of Lemma [3.5l For e > 0, 
let 

H{t,x,Keim),p) 


(3.33) 

where 

It follows that 


H^{t,x,m,p) = 


1 + eH{t, X, Ke(m),p)h“^)/^ 
1 


K^{m) := min < max{m, e}, 


(3.34) 


VpH^{t,x,m,p) = 


(1 + ^H{t,x,Ke{m),p)^'^ ^'>/'')VpH{t,x,Ke{m),p) 


(1 + eH{t, x, Ke(m),p)(''“^)/^)2 

so that by (j2.2j) and ()2.3j) we have 

\VpH{t,x,K^{m),p)\ 


\VpH^{t,x,m,p)\ < 


C. 


/ / X o/ <—{Ke{m) ^ + Ke{m) )<C{e). 


We also set := min{/, d}, := minl^r, d}. Note that Hf:{t,x,m,p), f^{t,x,m), and gf:{x,m) 

satisfy the hypotheses of Lemma 13.51 

Additionally, we claim also satisfies a form of the estimate (I3.30|) . which together with Lemma 
[Q is sufficient to give the a priori bounds from Corollary 13.31 From (j2.4p we see that 

(r - l)H{t,x,Ke{m),p) 


(3.35) 


VpH^{t, x,m,p) ■ p — H^{t, x, m,p) > 


(1 + eH{t, X, «;e(m),p)h“i)A)2 


-c. 


while from (j2.2j) and (j2.3j) we get 

(3.36) ^ 

< C{K^{m)^ + Ke(m)“'^)^^/^^“^^(1 + H^{t, X, m,p)) 

using the fact that r/(r — 1) > 2. We deduce 

(3.37) 

\VpHf:{t, X, + VpH^{t, x,m,p) ■ p — H^{t, x, m,p)), 

which gives an estimate of the form (13.3Up with an adjusted value of A. (Such an adjustment of the 
parameter A does not change the a priori estimates obtained from Corollary 13.31 1 
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Now let {u^,me) be bounded solutions of with H, /, g replaced by respectively. We 

claim that, up to a subsequence, {ue,m^) converges to a weak solution {u,m) of dni). Recall that 
we are assuming T is small enough so that Corollary 13.31 is valid. 

From Corollary 13.31 we see that —dtUe — Au^ is bounded in L^{Q) and that Ue{T) is bounded in 
L^(f7). This implies that and are relatively compact in L^{Q), so for some u G L^{Q) we 
have, up to a subsequence, ^ u strongly in L^{Q), —)• Vue strongly in L^{Q) as well as 

pointwise almost everywhere. 

On the other hand, we have that satisfies (13.4p . (13.5p . and (13.61) uniformly in e, that is, 

||"ie||oo + \\'>TT'e\\L‘^{0,T-,H^(n)) + l|9f”^e||L^(0,T;Vy-l’^(f2)) ^ C. 

From standard compactness results (cf. [15]) we have for some subsequence that strongly 
converges in L‘^{Q) and almost everywhere on a subsequence to some m. Moreover, by part 2 of 
Corollary 13.31 we have that \\VpHf:{t,x,me,'Vue)\\r/(r-i) < C, and since r/(r — 1) >2 we can use 
Holder’s inequality to see that me|VpFfe(t, x, m^, Viie)|^ is uniformly integrable. We conclude that 

X, rUe, Vtte) —)• y/mVpH{t, X, rn,'Vu) strongly in L^{Q). 

Applying [TH Theorem 6.1] we obtain —)• m in (^([O, T]; L^(H)) and that m is a weak solution 

of the corresponding Fokker-Planck equation. 

It remains to show that u solves the Hamilton-Jacobi equation. Since ^ m and Vttg —)■ Vrt 
almost everywhere, then by Fatou’s Lemma we have 

/ H{t,x,m,'Vu)dxdt / He{t,x,me,Vue)dxdt < C 

Jq "^0 Jq 

by (|3.26p . By hypothesis p2.2p together with the fact that m and m~^ are bounded, we see that 
\\Vu\\r <C SO that u G L''(0, T; Wb’'(H)). 

Now it can be shown using (I3.27P that ge{-,me(T)) and /e(-,-,me) are equi-integrable. Indeed, we 
have 

/ [ gL + j[ gei-,me{T))me{T) + C [ me(r). 

Je je ^ Je Je 

Taking the measure of E to zero and then L to infinity we see that the right-hand side goes to zero, 
as desired, so ge{-,rnt.{T)) is equi-integrable. The argument for is analogous. Given that 

me{T) —)■ m[T) and —>• m almost everywhere, it follows that 

5 e(-,me(r)) ^ g{-,m{T)) in L^(H), /e(-,-,m^) ^ f{-,-,rn) in L^{Q). 

As for the Hamiltonian term, because we have that and are uniformly bounded, H^{t, x, m^, V«e) 
therefore has uniformly subquadratic growth in by hypothesis p2.3p . We can now invoke pre¬ 
vious results on parabolic equations (see for instance m and references therein) to conclude that 
H^{t,x,m^,Vue) strongly converges in L^{Q) to H(t,x,m,Vu) and that u is a weak solution to 
the Hamilton-Jacobi equation. 

Finally, suppose we have that fi and gL are bounded for all L > 0. Then by p.jp it follows that 
f{m) and g{m{T)) are bounded. So by classical results on quasilinear parabolic equations, u is 
bounded and locally Holder continuous (see any text on parabolic equations, e.g. m)- 

This completes the proof of Theorem 12.21 

Remark 3.6. Naturally, we could show that the solution is even smooth, if we assume that H is 
smooth enough. 
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4. Uniqueness 


The proof of Theorem 12.41 requires two main ingredients. First, we prove that the estimates of 
Proposition 13.11 are valid for all weak solutions of the Fokker-Planck equation, and not just those 
with sufficient regularity. Second, we formulate the structure conditions on the Hamiltonian which 
allow us to prove uniqueness. 


4.1. Regularity of weak solutions. To obtain existence of solutions, it was sufficient to have a 
priori bounds on a sequence of smooth solutions, from which we showed that this sequence was 
compact. On the other hand, in order to establish uniqueness of weak solutions without restricting 
their definition, it is necessary to show that they are regular enough to perform integration by 
parts. For this we will prove an extension of Proposition 13.11 


First, we give a weaker definition of solutions of the Fokker-Planck equation (|3.ip - (l3.2l) . Suppose 
6g L2(Q). 

Definition 4.1 (Weak solutions of the Fokker-Planck equation). A function m € L^{Q)+ is a weak 
solution of (I3T])-(I32]) provided that 

m|6|2 gL^(Q) 

and 

/ — + / mo4>{0) 

J Q Jfl 

for every cf G C'“([0,r),H) satisfying the corresponding boundary conditions taken from ()2.8p . 


The above definition, taken from m, has two interesting properties. One is that, thanks to 
the condition that m|6|^ G L^{Q), one can show weak solutions defined in this way are unique. 
The second, which is of more interest to us, is that weak solutions are equivalent to renormalized 
solutions. 

Define, for k > 0, the truncation function 
(4.1) rfc(s) := min(A:, max(s, —fc)). 

A renormalized solution of (j3.1|) is defined as a function m G L^{Q)+ such that 

Tk{m) G L‘^iO,T;H^{Q)) 'ik > 0. 

lim ^ ff |Vmp = 0 

N^oo N JJN<\m\<2N 

f dtS{m) — AS{m) — div{S'{m)mb) = —5"(m)|Vmp — S"{m)mb ■ Vm 
\ 5'(m)(0) = S{mQ) 

for every S G VF^’°°(M) such that S' has compact support. Part (iii) is understood in the sense of 
distributions, again choosing test functions according to the specified boundary conditions. 

The equivalence between weak and renormalized solutions of (|3.ip will be crucial in the proof of 
the following: 

Proposition 4.2. Assume the same hypotheses as in Proposition 1,9.il except that m G L^{Q)^ is 
a weak solution of ()3.ip - ()3.2p . Then the conclusion of Proposition\3A\ holds. 


in Q 
in D 


(4.2) 


(i) 

(m) 


{Hi) 
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Proof. By [141 Theorem 3.6], m G (^([O, T]; L^(r2)) and is a renormalized solution. Let e > 0 be 
small, > 0 be large. Consider the auxiliary function 


(4.3) S{m) := < 


which satisfies S 
that S" > 0 we get 


e if m < 0 

(m + e)~‘^ if 0 < m < N 

{N + e)-'' + f (iV + e)-‘>-^[{m - iV - l)^ - 1] if iV < m < + 1 

[n + e)-^ - ^{N + ifN + l<m 

and S' has compact support. Now integrating (|4.2j) and using the fact 


(4.4) f S{m{t))P- f f S'"(m)|Vmp < f S'(mo) + - f [ S"{m)m'^\b\‘^. 

Jn 2 Jq Jq Jq 2 Jq Jq 

This implies 


(4.5) - f {m{t) + e) ^ ^ ^ ff (m + e) ^jVmp 

Q Jm{t)<N 2 J Jm<N 

<- [ m„-’+ |Si|(JV+ £)-« +// (m+ e)-»-2“|6|" + LjV+ £)-’-'(« +1)2 /■ |(,|2. 

QJn ^ JJm<N ^ Jq 

We may assume q is large. Let N —)• oo to obtain 

(4.6) - f {m{t) + e)~^ + ^ ^ ^ f [ (m + e)“'^“^|Vmp < - f ^ ^ ^ f /" (m + e)“^|6p. 

Q Jn J Jo Jn ^ Jn 2 Jq Jq 

Then we may follow the proof of Proposition 13.11 to see that (j3.5p holds. The proof of (|3.4p is 
analogous (replace —q with q). 


Finally, to see that (|3.6I1 holds, note that we have 0 < m < N almost everywhere for some N > 0. 
Define a new auxiliary function S G 1F^’°°(M) such that S' has compact support, and such that 
S{r) = for r G [0, A^j. Then (14.4p implies 


(4.7) 


^ f m{tf + \ [[ |Vmp f S{mo) 

J J J J 


+ 


m 



0 Jn 


\bf 


We conclude that m G L^(0, T; L7^(D)), from which (13.61) follows. 


□ 


Remark 4.3. By Definition \2.1l the result we have just proved implies that the conclusion of 
Proposition \3.1\ applies to m for any weak solution {u, m) of dni). 


Remark 4.4. One may consider proving Proposition \4.2\ by approximating weak solutions of the 
Fokker-Planck equation with smooth ones, appealing to Provosition \3.1\ and then passing to the limit. 
The difficulty one then encounters is to show that one can find such a smooth approximation such 
that (ITSI) also holds. For instance, if one takes an arbitrary smooth approximation b^ of the vector 
field b, then one risks losing the estimate p.3p . It turns out that using the theory of renormalized 
solutions is a much more tractable way of extending the essential estimates to weak solutions. 


4.2. Proof of Theorem 12.41 Now that we know weak solutions have some extra regularity, it will 
be possible to prove uniqueness using integration by parts. Before proceeding, we need a technical 
lemma. Recall that Assumption (12.9j) is in place, which allows us to prove the following: 

Lemma 4.5. Assume that (12.9|) holds. Then for all mi, mo > 0 and allpi,po G we have 
(4.8) {miVpH{mi,pi) - moVpH{mo,po)) ■ {pi - po) - {H{mi,pi) - H{mo,po)){mi - mo) > 0 
where strict inequality holds if pi po- 
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Proof. For 9 G [0,1], set pe = Po + 0{pi — po) and me = mo + 9{mi — mo), and let 

(j)[9) = {mgVpH{m 0 ,pg) - moVpH{mo,po)) ■ {pi - po) - {H{mg,pg) - H{mQ,pQ)){mi - mo). 
We compute 

f)\e) = -dmH{mg,pg){mi - mof + meVlH{me,pg){pi -po,Pi - Po) 

+ mg{mi - mo){pi - po) ■ dm'^pH{mg,pg). 

Use (12.9p to see that 4>'{9) > 0 for all 9 G [0,1]. Since (/)(0) = 0, it follows that (j){l) > 0, which 
is what we wanted to prove. In the case where pi po, it follows that pg = 0 for at most one 
9 G [0,1], so the above calculation shows 4>'{9) > 0 for all but at most one value of 0 G [0,1]. We 
conclude that (f{l) > 0, as desired. □ 


Now we prove Theorem 12.41 


Let (tti,mi) and (tt 2 ,m 2 ) be two solutions. Introduce the following sequence of auxiliary functions: 


(4.9) 


Sn{r) = nSi , Si{r) = J S[{t)dt, 5((r) 


1 if jsj < 1, 

2 — jsj if 1 < jsj < 2, 
0 if |s| > 2. 


Then for i = 1,2 and for /c > 0 we have that Sk{ui), Sk{u 2 ) G L^(0,T; Fl^(fl)) n L°°{Q), and Ui 
satisfies the renormalized equation 


(4.10) 


-dtSkiui) 


ASk{ui) + S{,{ui)H{t,x,mi,Vui) 




Sk{ui){T) = Sk{g{mi{T))). 


in the sense of distributions. Since by Proposition 14.21 we have mi,m 2 G f\L°°{Q), 

we can use mi — m 2 as a test function in (I4.10p to get 


(4.11) 

- [ [5^'(ui)|Vui|2 
Jq 


5’fc(M2)|Vu2p][mi-m2]+ [ [Skig{mi{T)))-Sk{gim 2 {T)))][mi{T)-m 2 iT)] 

Jn 


+ [ [S'kiui)f{mi)-Sk{u2)fim2)][mi-m2] 

JQ 

[ [Sk{ui)H{mi,Vui) - Sk{u2)H{m2,Vu2)][mi - m 2 ] 

Q 

+ [ [miVpiL(mi, Vtti) - m2VpiL(m2, Vn2)] • V[5fc(«i) - 5fc(«2)] = 0. 
JQ 


As A: —)■ 00 , S'f,{ui)\Vui\'^ —)• 0 strongly in L^{Q), for z = 1, 2. Since mi and m 2 are bounded, 
- [ [Sk{ui)\Vui\‘^ - Sk{u2)\Vu2\‘^][mi - m 2 ] ^ 0, k 

Jo 


00 . 


As for the other terms, from the dehnition of weak solution we know that, for z = 1, 2, LA(mj, Vztj) G 
L^{Q) which implies VpLI(mi, Vzzj) G by the assumptions on H and the fact that m* 

and m~^ are bounded. On the other hand, Vzzj G L'"{Q). So using the fact that S'(, —)• 1 and 
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Sk{u) —)• u and letting A; —)• oo in (I4.1ip . we get 

(4-12) [ [/(mi) -/(m 2 )][mi - m 2 ] + [ [g{mi{T)) - g{m 2 {T))][mi{T) - m 2 {T)] 

J Q Jfl 

+ / [miVpi/(mi, Vni)-m2VpiA(m2, Vn2)]-V[tti-n2]- / [-fA(mi, Vui)-//(m2, VM2)][mi-m2] = 0, 

JQ JQ 

Then we use Lemma 14.51 and the fact that / and g are strictly increasing to deduce that mi = m 2 , 
Vtti = Vu 2 , and ui{T) = U 2 {T), from which it follows that ui = U 2 as well. This completes the 
proof. 


5. Mean field type control 

Consider now the following system, which can be thought of as a companion to dni). 

(i) —dtu — Au + H(t,x,m,Vu)+mdmH{t,x,m,'Vu) = f{t,x,m) inQ 

(5.1) (m) dtm — Am — div {mVpH(t,x,m,Vu)) = 0 inQ 

{in) u{T,x) = g{x,m{T,x)), m(0,x) = mo(x) in 

The difference between (15.ip and (|l.ll) is the addition of the term mdmH{t, x, m, Vn), which makes 
the system into a condition of optimality for a control problem. See, for instance, [T]. Weak 
solutions of (15.11) can be defined in a way analogous to Definition 12.II for mean field games. 

Here we describe the mean field type optimal control problem solved by solutions to (|5.1I) . Let the 
Lagrangian L = L{t, x, m, b) be given by 

L{t, X, m, b) = supp ■ b — H{t,x, m, —p) 
p 

Let K, be the set of all pairs (m, w) G (L^(0, T; Lf^(D))+ n L°^{Q)) x such that 

dtm — Am + div(u;) = 0, m(0) = mo 

holds in the sense of distributions. Define 

L(t, X, m, w) = mL (t, x, m, — ) , {t,x) ^ Q, m > 0, tc G M'^. 

V m/ 

Observe that w ^ L{t, x, m, w) is obtained by taking the Legendre transform of p 1 —)• mH{t, x, m, —p). 
In this section we consider the optimal control problem of minimizing the objective functional 

(5.2) J{m,w) = / / L{t,x,m{t,x),w{t,x)) + F{t,x,m{t,x))dxdt + / G{x,m{T,x))dx 

Jo Jn Jn 

over all (m, w) G JC. Here F and G are given by 

pm pm 

F{t,x,m) = / f{t,x,s)ds, G{x,m) = / g{x,s)ds. 

Jo Jo 

The assumptions of Section 12.11 are still in force. We also add the following hypotheses: 

(1) We will assume in addition to the previous assumptions that for all (t,x) G Q, the map 
{m,p) i-A F[{t,x,m,p) is on (0, 00 ) x M'^. 

(2) We assume also that for all {t, x) G Q, the map (m, w) i-A L{t, x, m, w) is G^ and convex on 
(0, 00 ) X and that (|2.2p holds for dmH in place of H. 
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(3) Finally, we add the assumption that either {m,w) f-)- L{t,x,m,w) is strictly convex or else 
w !-)■ L{t, X, m, w) and m e-)- F{t, x, m) are both strictly convex. 

Let us remark that this last assumption on strict convexity is required for uniqueness of solutions. 
Existence of minimizers will hold even if the convexity is not strict. 

Remark 5.1. An example of a Hamiltonian H satisfying the above assumptions is the canonical 
example 

H{t,x,m,p) = ^ 

whenever a > 0 is any non-negative number. In particular, we remark that uniqueness holds for 
the mean field type control system 15.11 for a much larger class of Hamiltonians than for the mean 
field game (HH). 

In addition, we make the following deductions concerning F and G: 

• Since m >->■ f{t,x,m) and m >->■ g{x,m) are continuous and nondecreasing, it follows that 
m e-)■ F{t, X, m) and m i—>• G{x, m) are and convex. Moreover, it follows from (12.ip that 

(5.3) Fx,(t,x):= sup F{t,x,m), Gl{x) := sup G{x,m) 

m£[0,L] mg[0,L] 

are both integrable for all L > 0. 

• In addition, we remark that since / and g are bounded below, we have 

(5.4) F(t,x,m) > —Gm, G{x,m) > —Cm 
for some constant C > 0. 

Because L is strictly convex in re, we note that the equivalence of the hrst-order optimality criteria 

( 5 - 5 ) 

— p ='VwL{t,x,m,w) L{t,x,m,w) + mH{t,x,m,p) =—p ■ w —w = m'VpH{t,x,m,p). 

In particular, if any of these criteria hold, then by taking partial derivatives in m one obtains 
H{t,x,m,p) + mdmH{t,x,m,p) = —dmL{t,x,m,w). 

We now come to the main result of this section. 

Theorem 5.2. Under the assumptions above, and given T small enough, there exists a unique 
solution {u, m) of dSI]). The pair [m,—VpH{t,x,m,Vu)) is the unique minimizer of J in 1C. 

The proof of Theorem 15.21 comes in two Lemmas. 

Lemma 5.3. There exists a unique minimizer in 1C of J{m,w). 

Proof. Let {mn,Wn) be a sequence in 1C such that J{mn,Wn) converges to mi(m,w)£!C J- Set bn = 
Wnimn- Note that from (12.3p and the definition of L we can deduce that 

(5.6) L{t,x,m,b) > ' ' — ^-C{m^ + m^). 

Since {mn,mnbn) is a minimizing sequence of J, then using the estimates on F and G it follows 
that 

I mnL{t,x,mn,bn) < c -I C mn + C mn{T) < C\\mo\\i. 

Q J Q J Q. 
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(5.7) 




From (15.6p we deduce that, for some ?? > 0, 

(5.8) \bn\ < C{mPn + + Vn)-, := mnL{t, x, nin, bn), 

where < C by (j5.7p . 

Appealing to Proposition l3.ll it follows that mn is bounded in L^(0, T ; (fl)) and that mn and 

are bounded in L°°. Passing to a subsequence, we have that nin converges weakly in L^(0, T ; 
and pointwise almost everywhere to a function m G L°° such that 1/m is also in L°°. Moreover, 
(j5.6p and (j5.7p now imply that bn is a bounded sequence in Passing to a subsequence we 

conclude that bn converges weakly in to some b. 


We claim J{m,mb) < liminf^^oo J{'nin,'ninbn)- First, notice that for C > 0 large enough, 

/ F(t, x, m) + Cm < liminf / F{t,x,mn) + Crun 
JQ JQ 

by Fatou’s Lemma. By taking a subsequence we have m„ —)• m strongly in L^, so we get 

(5.9) / F(t,X,m) < liminf / F{t,x,mn)- 

Jq Jq 

A similar argument shows that 

(5.10) / G(x,m) < liminf / G{x,mn)- 

Jn Jq 

Next, because L is convex in (m,w), 


(5.11) 


/ Z(t, X, m, m6) < liminf / L{t,x,mn,mnbn)- 
Jq Jq 


Putting (15.9p . p5.10p . and ^5.111) together we see that J{m,mb) < liminf„_,.oo J(mn, rc^), which 
implies that J{m,mb) is the minimum over JC of J. 


To see that (m, mb) is unique, it suffices to note that by the assumptions given above, J is strictly 
convex. □ 


Lemma 5.4. If (m, w) G 1C is the minimizer of J, then [u, m) is a solution of (EH), where 
u G C([0, T]; L^(n)) solves the adjoint equation 

(5.12) - dtu - Au = f{m) + dmL{t,x,m,w), u{T) = g{m{T)). 

Conversely, if {u, m) is a solution of EH, then {m,—mVpH{t,x,m,Vu)) is in K, and is the 
minimizer of J. 


Proof. Let (m, w) G JC he the minimizer of J. By the proof of Lemma [5.3l m and 1/m are bounded. 
It follows that J is Gateaux differentiable at {m,w). To see this, observe that for any (mi,u>i) G JC 
we have that (m,w) = (mi —m,wi — w) is a pair in L^{0,T; H^{Cl)) x L^/(^“^)(Q) satisfying 

(5.13) dtfh — Am + div(t(;) = 0, m(0) = 0. 

We calculate 


(5.14) Wm.-Jim + 6m,w + 6w) 

e^o 9 


/ / dmL{t,x,m,w)m+ VwL{t,x,m,w) ■ w + f{t,x,m)m + / g{x,m{T))m{T). 

Jo Jn Jn 
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The right-hand side is well-defined: using (12.2p . (I5.6p . and (15.5p together with the fact that m and 
1/m are bounded we arrive at 

\VwL{t, X, m,w)\ < + C, 

which implies VwL{t,x,m,w) G because w G on the other hand, fh G L°°{Q) by 

Proposition 14.21 dmL{t,x,m,w) G L^{Q) by (|2.2p applied to dmH, and f{t,x,m) and g{x,m{T)) 
are integrable because m is bounded. 

Introduce the adjoint variable u G (^([OjT]; L^(fl)) by taking the solution of p5.12p . where we note 
that the right-hand side is in L^{Q). Then we have, by (15.131) . 


(5.15) 

\\m-J{m + 9m,w + 9w) = / / {—dtu — Au)m + VujL{t,x,m,w)-w+ / g{x,m{T))rh{T) 

^^0 9 Jo Jn Jq 

= / u{dtm — Am)+ VujL{t,x,m,w) ■ w = / / iVu + V^L{t,x,m,w)j-w. 

Jo Jn Jo Jn^ ' 

By optimality of (m, w), the value of this Gateaux derivative is zero. Since w is arbitrary, it follows 
that 

—Vu = VwL{t, X, m, w), 

and so by (15.511 we get 

(5.16) H(t, X, m, Vtt) -|- mdmH{t, x, m, Vu) = —dmL{t, x, m, w). 

Therefore by substitution u solves 

—dtu — Au + H{t, X, m, Vu) -|- mdmH(t, x, m, Vu) = f(t, x, m), u{T) = g{x, m{T)), 
while m solves 

dtm — Am — div{mVpH{t, x, m,p)) = 0, m(0) = mo, 

as desired. 


Finally, we show that solutions of (15.ip are equivalent to minimizers of J. Let (u, m) be a weak 
solution of m- Set w = —m'VpH{t,x,m,p). Note that {m,w) G JC. Following the reasoning 
in the previous step, we obtain that Vu = —VwL{t,x,m,w) and that (|5.16l) holds. Then by the 
calculations in (15.141) and (|5.15l) we get 


(5.17) 


lim - Jim + 9m, w + 
e^o 9 


9w) = f j ( 
Jo Jn ^ 


Vu -|- VyjL{t, x,m,w)j ■ w = 0 


hence the Gateaux derivative of J at (m, w) is zero. By the strict convexity of J, it follows that 
(m, w) is a minimizer of J. □ 


Proof of Theorem A5.Jl By Lemmas 15.31 and r5.41 there exists a solution (u, m) to (15.11) . On the other 
hand, if (u, m) is another solution to m, then (m, — VpLl(t, X, m, Vu)) is the unique minimizer 
of J, so m = m. Since u and u both satisfy the adjoint equation (I5.12p . it follows that u = u as 
well. □ 


To conclude, we remark on some possible generalizations of Theorem 15.21 The existence of min¬ 
imizers for the functional (j5.2p will hold under much broader assumptions than those presented 
above. In this article the hypotheses leading to the use of Proposition 13.11 ensure the Gateaux 
differentiability of (15.2p . allowing a classical optimal control argument yielding the solutions of 
(j5.1|) . For cases where the a priori bounds in Proposition 13.11 are not applicable, one may instead 
use the strategy adopted in [5] using Fenchel-Rockafellar duality. Accordingly, one seeks to prove 
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the existence of minimizers both for the optimal control problem min(^ „,)g^ J{m,w) and its dual, 
the latter of which is solved in a suitably relaxed setting. In this way one overcomes the lack of 
differentiability of J{m,w) due to its singular dependence on m. This method can also be used to 
construct solutions of (|5.ip under less restrictive conditions on the Hamiltonian. We can expect 
this to be the subject of future research. 
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